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^  INTRODUCTION 

Along  with  the  mean  velocity  that  takes  a  land  or  water  vehicle  from  one  place  to 
another,  there  are  six  types  of  oscillatory  motions,  which  are  rotational-roll,  pitch,  and  yaw; 
and  translational-surge,  sway,  and  heave.  These  six  types  of  motion,  as  well  as  the  mean 
vehicle  velocity,  cause  linear  velocities  that  add  to  the  muzzle  velocities  of  projectiles  fired 
from  guns  located  on  a  moving  vehicle.  ~~ 

Jr 

The  stabilization  system  of  a  gun  can  compensate  for  part  or  all  of  the  angular  devi¬ 
ations  of  the  aimpoint  caused  by  roll,  pitch,  and  yaw.  However,  the  linear  velocities  imparted 
by  these  rotations  as  well  as  the  translations  are  not  accounted  for  by  the  stable  element. 
Moreover,  there  are  some  fire  control  systems  in  use  that  do  not  account  for  the  same,  and 
so  in  these  systems,  the  unaccounted  velocities  affect  gun  fire  accuracy.  This  report  provides 
a  method  to  model  the  errors  caused  by  these  linear  velocities,  if  uncompensated  by  the  fire 
control.  Given  some  assumptions,  the  formulas  derived  are  exact  with  absolutely  no  approxi¬ 
mations,  and  they  can  be  utilized  by  writing  a  computer  program. 

This  work  was  necessary  since  the  operational  requirements  of  the  LVT(X)  amphibian 
vehicle  included  the  firing  of  its  weapon  from  water  near  the  shore.  Also,  the  Joint  Munitions 
Effectiveness  Manual  (JMEM)  Surface-to-Surface  Delivery  Accuracy  Working  Group  needed  a 
method  to  model  gun  fire  errors  due  to  firing  from  a  moving  platform. 


METHODOLOGY 

In  order  to  derive  the  model,  equations  that  calculate  the  linear  velocity  in  the  (North, 
West,  up)  frame  due  to  rotational  motion  are  explained.  Transformations  are  performed  on  all 
other  motions  of  the  vehicle  and  muzzle  velocities  to  express  them  in  the  above  frame.  The 
velocities  are  then  added  to  the  muzzle  velocity,  and  the  azimuth  and  elevation  formulas  for 
the  new  velocity  are  derived  from  the  sum.  The  errors  are  calculated  by  comparing  the  new 
azimuth  and  elevation  with  the  initial  ones  from  muzzle  velocity.  A  computer  program  using 
the  methodology  of  this  report  is  given  in  Appendix  A. 


LINEAR  VELOCITY  DUE  TO  ROTATIONAL  MOTION 

It  is  assumed  that  the  gun  barrel  is  pivoted  on  the  trunnion.  The  stabilization  system 
can  adjust  the  gun  barrel  to  compensate  for  angular  deviations  caused  by  vehicle  rotational 
oscillations,  but  it  cannot  adjust  for  the  linear  velocities.  Since  the  pivot  of  the  gun  barrel  is 
assumed  to  be  at  the  trunnion,  the  linear  velocity  imparted  on  the  projectile  from  vehicle 
motion  is  equal  to  the  linear  velocity  on  the  trunnion.  This  is  because  the  gun  barrel  and 
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muzzle  do  not  rotate  with  the  vehicle  (due  to  the  stable  element)  while  the  trunnion  does.  A 
small  additional  velocity  imparted  while  the  projectile  travels  down  the  gun  barrel  will  be 
discussed  in  Appendix  B. 

Assuming  that  the  centers  of  rotation  of  roll,  pitch,  and  yaw  coincide,  define  X'  pointing 
towards  the  fore,  Z'  pointing  up,  and  Y'  pointing  in  the  direction  that  forms  a  right-handed 
system.  Also,  X'  is  in  the  centerline  of  the  vehicle,  and  the  point  (0,0,0)  is  in  the  center  of 
rotation. 


An  equation  will  now  be  derived  for  the  linear  velocity  with  respect  to  a  (North,  West, 
up)  frame  imparted  by  the  rotational  motion. 

Let  r  =  roll,  p  =  pitch,  h  =  heading;  assuming  the  angle  of  the  mean  vehicle  velocity 
vector  is  constant, 

dh  /d  YAW  \ 

Also,  let  A  =  roll,  B  =  pitch,  C  =  heading  coordinate  rotations.  The  positive  direction  of  roll 
and  heading  are  opposite  of  the  conventions  in  a  right-handed  system.  The  'Vs  in  the  following 
frames  show  that  the  axis  is  not  in  the  (North  West,  up)  frame. 

/’  0  0  \ 

A  =  I  0  cos  r  sin  r  I  roll 

\  0  -sin  r  cos  r  / 
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pitch 


Multiply  CBA  by  the  position  vector  (R')  of  the  trunnion  in  the  primed  (platform)  frame 
gives  its  position  vector  in  the  (X,  Y,  Z)  *  (North,  West,  up)  =  (N.W.U)  frame,  which  is 
assumed  to  be  fixed  (Equation  (2)).  Hence, 


where  is  the  roll,  R^  is  the  pitch,  and  R^  is  the  heading  and  yaw  axes. 

The  linear  velocity  with  respect  to  (X,Y,Z)  is  simply  the  derivative  of  R  =  R(t)  as  shown 
in  Equation  (3). 
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In  order  to  calculate  the  derivative  on  the  right-hand  side  of  Equation  (3),  the  following 
definition  must  be  made  by  Equation  (4).  Given  matrix 

M  «  (My)  <4> 

Equation  (5)  is  the  definition  of  the  derivative  of  a  matrix,  which  is  the  derivative  of  every 
element  in  the  matrix. 


dM 

dt 


(5) 


Equation  (6)  is  a  mathematical  theorem.  Given  matrices  A,,  A2,  ...  An  and  vector  R 


dA, 


dR 


The  proof  is  given  in  Appendix  C. 


(6) 


From  Equations  (4)  and  (5),  (V^,  V  ,  Vf)  can  be  calculated  as  follows. 


dC  ,  dB  ,  dA  A  dR' 

=  —  BAR  +  C  —  AR'  +  CB  —  R  +  C  B  A  — 
dt  dt  dt  dt 


(7) 


Note  that  the  last  term  on  the  right  is  zero,  since  R'  is  assumed  fixed  on  the  vehicle 
coordinate  system.  Hence, 
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*  —BAR'  +  C  —  AR'  +  CB  —  R' 
dt  dt  dt 


where 

V  represents  the  linear  velocity  imparted  by  rotational  motion. 

The  matrix  derivatives  shown  are  obtained  by  differentiating  each  element  of  the  original 
matrices. 


f  1  0  0  \  /o 

-  -  -f  \( 

dt  dtj  0  cost  sin  r  |=|  0 

\  0  -sin  r  cos  r/  \  0 


dr  dr 

-  sin  r  —  cos  r — 
dt  dt 


-  cos  r 


/ .  tip 
/-smp- 

0 

cos  p 

dP\ 

/  dh 

(  dt 

,  dh 
cos(h)  — 
dt 

1  ° 

0 

0 

) 

dC  I 
dt  ~  1 

x  dh 
-cos (h)  — 
dt 

•  dh 
-sin  (h)  — 
dt 

\  dp 

X-cosp  - 

0 

-sin  p 

dp_  J 

dt/ 

V  ° 

0 

/ 

V 

An  alternate  equation  can  be  derived  to  calculate  the  linear  velocity  due  to  rotational 
motion  by  noting  that 

dR 

V,„,  -  -  =  W  x  R  (9) 

where  R  and  W  are  both  expressed  with  respect  to  a  fixed  frame  that  is  assumed  to  be 
(N,W,U).  From  previous  explanations,  it  is  known  that 


=  CBAR' 
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where  R'  is  the  same  as  the  one  that  was  already  defined.  W  must  be  found,  which  is  the 
rotation  with  respect  to  the  fixed  frame.  Roll,  pitch,  and  heading  are  given.  Only  heading, 
whose  negative  direction  is  up,  can  be  considered  in  the  fixed  frame.  The  roll  axis  can  be 
transformed  into  the  fixed  frame  via  rotations  with  magnitudes  equal  to  pitch  and  heading, 
respectively,  but  in  opposite  directions.  The  pitch  axis  is  parallel  to  the  earth’s  horizontal  axes 
and  at  an  angle  of  the  heading’s  magnitude  clockwise  from  the  West.  Hence,  the  pitch  axis 
can  be  transformed  to  the  fixed  frame  by  a  rotation  of  magnitude  equal  to  heading  but  in 
the  opposite  direction.  Hence, 


W  =  CB 


(11) 


LINEAR  VELOCITY  DUE  TO  TRANSLATIONAL  MOTION 

In  order  to  describe  the  linear  velocity  caused  by  surge,  sway,  and  heave,  a  precise  mathe¬ 
matical  definition  ot  these  motions  must  be  made.  Since  there  seems  to  be  some  ambiguity 
on  an  accepted  definition,  they  are  defined  in  the  next  paragraph. 

The  surge,  sway,  and  heave  are  oscillatory  translational  motions  defined  on  a  coordinate 
system  whose  origin  is  the  mean  location  of  the  center  of  rotation  that  is  assumed  to  be 
at  or  very  near  the  mean  location  of  the  center  of  mass.  The  X-axis  is  the  surge  axis  and 
points  toward  the  direction  of  the  mean  velocity  vector,  if  the  vehicle  is  moving  forward;  and 
negative  of  this  vector,  if  the  vehicle  is  going  backwards.  The  Z-axis  is  the  heave  axis  and  is 
pointed  up  with  respect  to  the  earth.  The  sway  axis  is  horizontal  and  points  in  a  direction 
that  completes  a  right-handed  system.  As  above,  assume  the  vehicle  mean  velocity  vector  is 
constant.  If  the  vehicle  is  not  traveling,  such  as  a  ship  floating  in  water,  the  X-axis  is  defined 
as  the  mean  location  of  the  centerline  of  the  vehicle. 

Given  the  translational  oscillations  in  the  above  coordinate  system,  these  motions  can  be 
expressed  in  the  (N,W,U)  frame  by  the  following  transformation: 


where  0  is  the  vehicle  mean  velocity  angle  from  the  North.  As  before,  let  positive  0  be  the 
opposite  of  the  convention,  that  is,  0  is  defined  as  clockwise  if  0  >0. 
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X  =  SURGE  AXIS  *  VELOCITY  VECTOR  AXIS  FOR  FORWARO  MOTION 


The  remaining  translational  motion  is  the  mean  vehicle  velocity  itself,  and  it  is  not 
oscillatory.  Given  the  mean  speed  and  $  as  defined  above,  this  velocity  in  the  (N,W,U)  frame 
is  as  follows: 


V 


mein 


vel 


I  Vm..n  v.l  I  (C0S«-  '  Sin*’  °> 


(13) 


ACTUAL  PROJECTILE  VELOCITY 

The  actual  projectile  velocity  is  the  vector  sum  of  the  given  velocity,  which  is  the  muzzle 
velocity  VmulJlt  plus  all  of  the  above  linear  velocities  imparted  to  the  projectile  by  vehicle 
motion  at  t  =  to  when  the  gun  is  fired.  In  the  (N,W,U)  frame 
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*  J 


m  unit 


I  v»«..u  I  (co»(H#)  °os( A,*),  -  cos (Elo)  sin (Af*),  V  sin(Elo)) 


where 


Elo  *  initial  gun  elevation 


A#*  =  Az(initial)  +  Heading 


At  (initial)  »  gun  bearing  with  respect  to  forward  of  the  platform 

Adding  all  the  various  velocities  at  t  =  to>  the  actual  bullet  velocity  Vb  is  as  follows: 

v  =  V  +  IV  1 

*b  muttle  L  v«hicl«  J  «  -  t 

=  V  +IV+V  +V  1  (15) 

muttle  I  rot  otc  trintl  mein  vet  J  t  *  t  '  ' 

o 

where  b  denotes  bullet.  V  .  is  often  assumed  to  be  compensated  after  the  first  round, 

ITlfin  V®l  1 

and  so  it  may  need  to  be  removed  in  a  computer  model  after  the  first  round. 

The  correspondence  of  the  terms  of  Equation  (15)  to  that  derived  from  Coriolis’  of 
classical  mechanics  at  t  further  demonstrates  the  truth  of  the  former.  Let  R  be  the  t 

o  o 

bullet  position  in  (N,W,U).  Let  R  be  the  fo  bullet  position  in  the  vehicle  frame  assuming 
for  now  R'  -  trunnion  position  at  t  and  L  be  the  vehicle  frame's  position  in  (N.W.U). 


From  Coriolis’  equations  in  two  forms, 


+  W  x  R' 


d'R:  dM  , 

M  -  +  -  R 


d  d' 

where  —  and  —  are  derivatives  in  the  (N,W,U)  and  (X',  Y\  Z')  frames,  respectively,  M  = 
dt  dt 

rotation  matrix  from  (X'.  Y\  Z')  to  (N,W,U)  and  R'  is  R'  in  terms  of  (X',  Y',  Z').  Hence, 


dR  d'R'  ,  dL 

—  =  -  +  W  x  R  + -  and  also 

dt  dt  dt 


$ 


’  *  \*  V  M  .  ' •  »  u9*' 
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dR 

dt 


d'R* 


dM  dL 

+  - R'  +  — 

dt  T  dt 


(20) 


d'R'  d'R'v 

Letting  R'  be  in  tenns  of  (N,W,U),  — ■ —  of  Equation  (19)  and  M  of  Equation  (20)  both 

dt  dt 

equal  Vmutll#  with  respect  to  (X',  Y',  Z')  at  t0  but  in  terms  of  (N.W.U). 

W  x  R'  *  Wx  (trunnion  position)  and 


—  R'  *  -  x  (trunnion  position),, 

dt  v  dt  v 

both  give  V  .  at  t  with  respect  to  and  in  terms  of  (N.W.U).  Note  that  M  is  C  B  A  and 

^  f  O  t  O 

dM 

—  *  CBA  +  CBA  +  CBA  as  before; 
dt 


dL 

dt 


m«tn  v#! 


+  Vo.«U  (rant  and 


dR 

dt 


=  V 


b 


Hence,  Equations  (19)  and  (20)  contain  exactly  the  terms  on  the  right  side  of  Equa¬ 
tion  (15)  at  t  . 

O 

Since  the  gun  is  pivoted  on  the  trunnion,  the  above  holds  even  if  R'  ¥=  trunnion  position 
at  t0.  In  this  case,  one  adds  a  fixed  translation  to  (X',  Y',  Z')  so  that  the  new  origin  is  the 
center  of  rotation  of  the  bullet’s  initial  position.  The  vector  from  this  new  location  to  the 
bullet  position  is  still  R'  or  R^  and  so  all  terms  containing  R'  or  R^  are  unchanged.  Letting 

a  A 

L  *  new  L,  L  *  L  +  fixed  translation  implies  that 

dL  ^  dL 
dt  dt 

Letting  R  *  new  R,  6  *  L  +  R'  since  R  =  L  +  R',  and  so 

dR_dL+dfV  MdL  +  dfV  _  dR 
dt  dt  dt  dt  dt  dt 

Therefore,  assuming  R'  ¥>  trunnion  position  at  t  changes  nothing. 
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CALCULATION  OF  ERRORS 


El  =  sin" 1  — 

new  y 


(22) 


AZ* 

new 


*  COS' 


\Zv2  +  vJ 


if  Vy  <  0 


(23) 


AZn.w  =  -COS'' 


v/v7  TW 


(24) 


if  V  >0,  since  clockwise  is  defined  positive  in  this  application. 


See  the  following  diagram. 


v  =  w 


y  =  w 


—  X  -  N 


AZ* 

new 


=  AZ  +  Heading 

new  ^ 


(25) 


The  azimuth  and  elevation  errors  are  as  follows: 

e .  =  AZ*  -  AZ*  (26) 

A  L  o 
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€ 


El 


El 


El 


o 


The  distance  errors  for  a  flat  trajectory  are  as  follows: 


e\z  °  (RANGE)  <«AZ) 


(27) 


(28) 


e'E1  *  (RANGE)  (eE1)  (29) 

Hence,  we  have  modeled  the  gun  shot  errors  caused  by  linear  velocities  imparted  by  vehicle 
motion.  If  one  desires  a  more  precise  assumption  than  a  flat  trajectory,  other  formulas  not 
given  here  can  be  used  to  find  (e'AZ,  e'EI)  when  (eA2,  eEI)  as  calculated  above  is  given. 


MODELS  OF  VEHICLE  MOTION 

A  time  series  model  of  vehicle  motion  in  water  or  land  can  be  obtained  experimentally. 
The  following  is  a  simple  analytical  model  of  vehicle  motion  in  water,  which  is  implemented 
in  the  model  and  may  be  useful  in  some  applications. 


r  =  roll  *  MAX 


■  *(?') 

(? 

h  =  heading  »  MAX^  sin 


p  ■  pitch  =  MAXp  sin  ( —  t 


t  +  <P 


(30) 


(31) 


(32) 


where  T,s  are  periods  and  MAX,s  are  maximum  amplitudes,  <f>  »  angle  of  ship  velocity  vector 
(assumed  constant),  and  t  =  time. 


*  ■  ip  MAX  cos /^tj 


(33) 
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*■?  MAx-w,(xr) 


2w 


h*  V  Maxy.w  «»  Ur 


(?•) 


(34) 

(35) 


For  surge,  sway,  and  heave,  only  tHeir  derivatives  are  used,  that  is,  their  velocities.  So, 


v.»,„  “  7-  MAX  (SURGE)  cos 


(?•) 

V,„„  -  ~  MAX  (SWAY)  C0i  t\ 

It) 


hMV«  -j- 


A 

^  MAX  (HEAVE)  cos  (  ^ 

S 


(36) 

(37) 

(38) 


where  T,s  are  the  periods  and  t  »  time. 


IMPLEMENTATION  OF  THE  MODEL 


A  computer  program,  AIMPT,  was  written  in  FORTRAN  that  gives  the  aimpoint  error 
(azimuth  and  elevation)  due  to  linear  velocities,  where  the  linear  velocity  imparted  by  rotational 
motion  was  modeled  by  the  matrices  method.  Also,  platform  motion  itself  was  simulated  in  a 
program  called  SHIP  that  uses  a  sinusoidal  model.  A  program  called  H1TPROB,  written  by  the 
U.S.  Army  Materiel  Systems  Analysis  Activity*  was  used  as  a  driver  to  the  linear  velocity 
model.  This  program  simulates  the  probability  of  hitting  a  square  or  rectangular  target  by 
projectiles  fired  from  a  25mm  chain  gun  mounted  on  a  Bradley  Fighting  Vehicle.  Listings 
of  AIMPT  and  SHIP  are  given  below,  and  the  program  was  run  one  time  without  linear 
velocities  and  one  time  with  Unear  velocities.  The  Usting  of  the  program  is  in  Appendix  A. 
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Input  parameters  to  AIMPT  and  SHIP  were  as  follows: 


Gun  position 

20°  AZ  20°  El 

Bullet  muzzle  velocity 

134S  m/sec 

Bearing  at  vehicle  from  north 

60* 

Maximum  roll 

0.0873  rad 

Maximum  pitch 

0.0873  rad 

Maximum  yaw 

0.0349  rad 

Trunnion  location  from  center 

of  rotation  (1.,  1.73205,  1.73205)  in  m  (meters) 

Roll  period 

2  sec 

Pitch  period 

2  sec 

Yaw  period 

5  sec 

Vehicle  forward  speed 

4  m/sec 

Maximum  surge  amplitude 

0.1  m 

The  output  is  as  follows 

Range  2000  m 

Target  size  2.286  x 


Maximum  sway  amplitude 

0.1  m 

Maximum  heave  amplitue 

0.5  m 

Surge  period 

10  sec 

Sway  period 

10  sec 

Heave  period 

2  sec 

Time  increment 

0.6  sec 

.286  ml 


Hit  probability  0.453 

without 

linear  velocities 


Hit  probability  0.348 

with 

linear  velocities 


As  seen,  the  probability  of  hit  has  been  noticeably  lowered  by  the  presence  of  linear  velocities. 


CONCLUSION 


The  model  explained  in  this  report  and  implemented  by  the  computer  program  described 
shows  that  linear  velocities,  imparted  by  gun  platform  motion  on  projectiles,  influence  the 
probability  of  hitting  a  target. 
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RECOMMENDATION 


In  modeling  gun  fire  accuracy,  consideration  should  be  given  to  the  effects  of  linear 
velocities,  imparted  by  platform  motion,  on  accuracy.  A  model  such  as  the  one  described  in 
this  report  or  actual  test  data  can  be  used  to  model  this  effect. 
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APPENDIX  A 
COMPUTER  PROGRAM 
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APPENDIX  B 

VELOCITY  IMPARTED  ON  A  PROJECTILE  AS 
IT  TRAVELS  THROUGH  A  GUN  TUBE 
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A  source  of  additional  velocity  not  yet  analyzed  in  this  report  is  the  velocity  imparted 
by  angular  motion  as  the  projectile  travels  through  the  gun  tube.  This  velocity,  V  ,  is  caused 
by  the  vehicle  angular  velocity  continuing  to  impart  linear  velocity  when  all  or  part  of  the 
moving  bullet  is  still  in  the  tube. 


Since  the  component  of  Vgunthat  is  parallel  to  Vb  does  not  totally  add  onto  the  pro¬ 
jectile,  only  the  component  of  Vgun  normal  to  Vfe  completely  adds  to  Vfe.  Call  this  component 
Vg.  Hence  |  Vgun|  =  max  |  Vg|  .  Since  the  stabilized  gun  tube  does  not  rotate  with  the 
vehicle,  the  magnitude  of  the  radius  of  rotation  on  the  tube  can  be  represented  by  Vrot 
defined  previously  in  this  report.  Hence,  given  that  t2  -  t,  represents  the  elapsed  time  from 
discharge  to  the  time  the  bullet  completely  exits  the  gun,  we  have  the  following  equation 


vrot(t)dt  =  vrot(t2)  -  vot  (t.) 


(A-l) 


This  velocity  is  generally  very  small  compared  to  V  ,  since  t2  -  t,  is  small  and  so 
V  .  (t2)  as  V  .  (t, ).  Hence,  neither  V  nor  V  is  entered  in  the  model. 

rot*7  rotvl/  g  gun 
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APPENDIX  C 

FORMULA  FOR  DIFFERENTIATING  A  PRODUCT  OF 
MATRICES  TIMES  A  VECTOR 
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d  E  n  dA.  n  dR 

Theorem;  —  n  A.  R  =  2  (A.  ...  -  ...  A  R)  +  II  A.  - 

dt  i=l  *  i=l  1  dt  n  i=l  1  dt 


Proof; 

Mathematical  induction  is  used. 


(B-l) 


(B-2) 


(B-3) 


(B-4) 


(B-5) 
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dt 


dRj 

+  a"J  "dt" 


AR  dA  dR 

—  -  —  R  +  A— 
dt  dt  dt 


-  r 


(B-6) 


(R-7) 


(case  for  n+1  matrices  if  true  for  n  matrices) 
Suppose 


d_ 

dt 


n 

n 

i«l 


n 


A.  R  •  I 
i-1 


n 

+  n  a. 

i“l 


dR 

dt 


(B-8) 


d  n+l 

7-  n 

dt 


A, 


R 


Without  loss  of  generality,  the  subscripts  can  be  renamed  as 
ie  jl,  ...  ,  n+lj 


(B-9) 

n  l  and  then  back  to 


£ 

dt 


n+1 

n  a.  r- 

i-1 


d\ 

m  — 

dt 

d\ 

m  ■  ■  ■ 

dt 

dA, 

n  • 

dt 


n  dA.  /  n 

n  Aj^r 

k»l 


dA  /  n  \  d  n 

— r  (  n  Ak  R )  +  ao  —  n  a.  r 

d<  \k-l  1  )  d‘  kM 

/  »  \  r »  dR  ;  /  dA 

\k-l  *  /  0  [_k-l  ^  d:  k-l\  dl 

[;,  A  *}  .:,(*•  A 


n 

n 

k-1 


n  AkR  +  Ao 


d \ 
dt 


n+1  /n+1  do  \  n+1  /  dA 

n  A.R  +(  n  A.  —  +  I  A  A  ...  —  ... 

i-2  '  V-l  ‘  d<  /  i-2  \  '  ‘  dt 


(B-10) 
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n+1 

dR 

dA, 

n+1 

n+1 

/  dA. 

n 

i°l 

A.  “7“ 
1  dt 

+  - 

dt 

n 

i«2 

AjR 

♦  2 
i-2 

(A»A»-  dt 

n+1 

dR 

n+1  / 

d^ 

\ 

n 

i-l 

A‘  dt 

+  2  ( 
i=l  \ 

dt 

V.R) 

l  true  for  n+1 

Thus,  by  mathematical  induction,  the  theorem  is  true. 
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